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A QUANTITATIVE  ANALYSIS  OF  ACOUSTIC 
CAVITATION  INCEPTION  IN  WATER 
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INTRODUCTION 

Many  attempts  have  been  made  to  explain  acoustic  cavitation 
inception  thresholds  but  have  achieved  little  success.  The  work  of 
Harvey,1  Barger,2  Strasberg, 3 and  Apfel4  is  well  known,  and  Holl,  5 
in  a recent  review,  lists  51  papers  concerning  the  inception  problem. 

A historical  survey  of  the  problem  will  not  be  given  here;  recent 
reviews5’6  are  quite  complete.  The  crux  of  the  problem  in  explaining 
cavitation  thresholds  is  the  proper  treatment  of  the  nuclei  that  induce 
cavitation.  Holl,5  for  example,  concludes  his  review  with  the  comment 
"It  is  very  apparent  that  the  understanding  of  the  onset  of  cavitation 
is  greatly  dependent  upon  knowing  the  form,  size  distribution  and  source 
of  cavitation  nuclei."  A renewed  interest  in  the  onset  of  cavitation  in 
liquids  has  been  generated  by  the  discovery  that  polymer  additives  that 
reduce  drag  also  reduce  cavitation.7'8  This  point  is  now  well  established 
and  is  of  considerable  interest  because  of  its  double  positive  effect. 

An  adequate  explanation  of  the  cavitation  reduction  problem  has  not  been 
given  but  attention  is  being  given  to  fluid  dynamics  rather  than  effects 
of  physical  parameters.  In  particular,  Ellis9  has  commented  that 
"an  explanation  of  the  reduced  cavitation  in  polymer  solutions  must  lie 
in  changed  flow  dynamics  rather  than  fluid  physical  properties". 

This  report  will  present  evidence  that  the  cavitation  inception 
threshold  can  be  explained  in  terms  of  nucleation  from  gas-filled 
cavities  in  motes,  that  is,  microscopic  particles  of  nonpolar  solids 
that  are  suspended  in  the  cavitating  liquid.  A theoretical  expression 
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for  the  cavitation  threshold  will  be  given  that  is  independent  of  the 
geometrical  properties  of  the  suspended  particle  and  that  correctly 
predicts  the  behavior  of  the  cavitation  threshold  • ' rh  temperature, 
frequency,  liquid  surface  tension  and  dissolved  gas  content. 


RESULTS 


A.  Theoretical  model 

The  model  that  will  be  used  is  one  originally  due  to  Harvey1  that 
has  been  modified  by  Strasberg3  and  more  recently  by  Apfel.1*  The  general 
approach  of  Apfel  will  be  followed;  however,  details  of  the  calculations 
will  be  given  here  for  completeness. 

It  is  assumed  that  there  exist  in  the  body  of  the  fluid 
numerous  particles  of  nonpolar  solids  that  contain  gas  that  is  trapped 
in  irregularities  in  the  particle's  surface.  For  simplicity  it  is  assumed 
that  the  irregularities  containing  the  entrained  gas  are  cones,  but  the 
approach  is  easily  generalized  to  wedge-shaped  cracks.  A model  mote 
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In  order  for  this  mote  to  be  an  effective  cavitation  nucleus,  it 
must  be  capable  of  stabilization,  i.e.,  it  must  not  allow  the  gas  to 
dissolve  in  the  normal  experimental  condition  of  a partially  degassed 
liquid.  This  will  occur  when  the  interface  is  concave  outwards  to  such 
an  extent  that  forces  on  the  interface  are  in  equilibrium,  viz., 

?h  - pv  + Pg  + -y  . 

where  Pjj  is  the  hydrostatic  pressure,  Pv  is  the  vapor  pressure  of  the 
liquid,  Pg  is  the  equilibrium  gas  content  of  the  liquid,  a is  the  surface 
tension  and  R is  the  radius  of  curvature  of  the  interface.  In  addition 
there  exist  requirements  on  the  advancing  contact  angle  otA  and  crevice 
angle  20; 

0 + | - aA  1 ti  . (2) 

It  is  appropriate  now  also  to  make  the  assumption  that  for  crevices  to 
be  effective  they  must  have  small  angles,  that  is,  cos  0 ~ 1.  It  is 
known10  that  nonpolar  solids  that  are  good  candidates  for  mote  material, 
e.g.,  parafin,  beeswax,  cellulose  acetate  have  advancing  contact  angles 
on  the  order  of  100°.  Thus,  for  Eq.  2 to  hold,  0 ~ 10°.  If  motes 
similar  to  these  exist  in  a degassed  liquid,  the  air-liquid  interface 
will  advance  until  the  above  equations  are  obeyed,  and  will  then  stop. 
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At  this  stabilization  position, 

a = R cos  (aA  - 3)  • (3) 

■I 

If  Eq.  1 is  solved  for  R,  the  half-width  of  the  crevice  can  be  expressed 
as 


p p— . p cos  (aA  - 3)  . (4) 

*v  *g 


If  we  assume  that  the  motes  are  all  nonpolar  solids  such  as  parafin, 
and  the  crevices  all  have  small  angles,  then  it  is  convenient  to  set  the 
contact  angle  relationship  equal  to  a constant  6.  Consider  next  a typical 
experimental  situation  in  which  the  liquid  is  subjected  to  a negative 
pressure  Pl-  The  interface  will  now  bow  out  and  take  the  shape  shown 
in  Fig.  2.  I^L 

->i  a k- 


Figure  2.  Nucleation  from  a mote. 
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The  force  balance  now  requires 


PL  = pv  + Ypg  - R.  • (5) 

•I 

The  factor  y has  been  added  to  account  for  the  fact  that  in  the  dynamic 
situation  the  diffusion  of  gas  across  the  interface  may  not  be  rapid 
enough  to  maintain  equilibrium.  Again,  R can  be  expressed  in  terms  of 
the  liquid-solid  contact  angle,  aR  as 

R'  = a/cos  (aR  - 6)  . (6) 


a has  already  been  set  by  the  initial  conditions  as  given  in  Eq.  4. 
The  radius  of  curvature  of  the  interface  is  then 

R'  = 2o 6_ 

PH  ~ pv  “ pg  cos  (aR  - 6)  ’ 

and  the  pressure  balance  on  the  interface  can  be  expressed  as 


P 


Li 


Pv  + Vpg 


(PH_“  pv  ~ pg)  cos  (aR  ~ $) 

6 


(8) 


If  the  negative  pressure  imposed  on  the  liquid  is  sufficient  to  cause  the 
interface  to  become  hemispherical  or  to  reach  the  receding  contact  angle, 
then  the  interface  will  advance  along  the  crevice  and  is  said  to  have 
nucleated.  Since  we  have  assumed  8 to  be  quite  small,  and  since  it  is 
known  that  aR  =»  40°  for  nonpolar  solids,10  the  requirement  for  nucleation 


8 


will  be  taken  as  the  condition  that  the  interface  achieve  the  receding 
contact  angle.  The  pressure  in  the  liquid  is  related  to  the  acoustic 
pressure  by  the  relation  = PL  ~ PH  » so  that  the  acoustic  pressure 
amplitude  required  for  cavitation  inception  is 


PA  =-(PH  ~ Pv  ~ YPg) 


(?H  ~ Py  - Pg)  cc  (CtR  ~ 3) 
6 


(9) 


This  is  essentially  equal  to  Eq.  3-13a  of  Ref.  4. 

In  order  to  apply  this  theory  to  an  experimental  situation,  it  is 
necessary  to  specify  the  contact  angle  relationships.  The  actual  procedure 
of  the  experimental  measurements  is  also  very  important.  To  test  the 
theory,  Eq.  9 will  be  applied  to  measurements  of  the  variation  of  the 
cavitation  threshold  with  gas  concentration,  temperature  and  surface 
tension.  The  variation  with  surface  tension  is  particularly  important  for 
it  requires  that  the  contact  angles  also  be  expressed  as  a function  of 
surface  tension.  In  the  experimental  procedure,  the  liquid  is  first 
degassed  to  a particular  value,  varying  amounts  of  surfactant  are  added, 
and  then  the  cavitation  threshold  is  determined.  The  initial  position 
of  the  interface  in  the  crevice  is  set  when  the  surface  tension  is  at  its 
initial  value,  and  when  surfactant  is  added,  does  not  change;  thus,  6 is 
not  a function  of  a.  As  the  surfactant  is  added,  the  surface  tension 
is  reduced,  but  Eq.  1,  which  specifies  the  conditions  for  the  equilibium 
of  the  interface,  shows  that  the  radius  of  curvature  will  then  reduce  to 
compensate  for  the  change  in  surface  tension.  The  receding  contact  angle 
relationship  will  be  a function  of  surface  tension,  however.  For  nonpolar 
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solids  such  as  parafin  there  is  a significant  hysteresis  between  the 
equilibium  contact  angle  and  the  receding  contact  angle.11  Hysteresis 
is  necessary  for  this  stabilization  mechanism  to  work,  of  course,  and 

must  be  accounted  for.  It  is  possible  to  write 

■I 

cos  (aR  - 3)  = cos  (otg  - ajj  - 3)  “ cos  [ae  - (oh  + 3)]  . (10) 

which  can  be  expanded  to  give 

cos  (aR  - 3)  = cos  ae  cos  (<}))  + sin  ae  sin  (<J>)  , (11) 


where  ae  is  the  equilibium  contact  angle,  ay  is  the  amount  of  hysteresis 
and  <}>  = ay  + 3*  Bargeman  and  Van  Voorst  Vader12  have  determined  that  the 
effect  of  surfactants  on  the  equilibium  contact  angle  can  be  expressed  by 


cos  ae  = -1  + C/a  , 


(12) 


where  C is  a constant  that  depends  on  the  surface  properties  of  the  solid. 
For  parafin,  C = 49,  for  beeswax,  C = 48;  we  have  chosen  C = 48.5.  The 
acoustic  cavitation  threshold  can  then  be  expressed  in  terms  of  measureable 
parameters  as 


pA 


(Py  - Pv  - YPg)  “ jp [cos  ((}>)  (C/a  - 1)  + 

sin  <f>  [l  - (§  - 


(13) 


10 


B . Application  of  theory  to  experimental  results 


The  incipient  threshold  was  measured  for  acoustic  cavitation  in 
distilled  water.  Various  hollow  piezoelectric  cylindrical  transducers 

were  used  that  were  open  at  one  end  and  closed  at  the  other  with  a 

■I 

tightly  stretched  stainless  steel  foil.  A typical  one  would  have  an 
inside  diameter  of  6.2  cm,  a height  of  7.5  cm  and  would  be  driven  in 
(r,0,z)  = (3,0,3)  normal  mode.  Higher  modes  restricted  the  cavitation 
to  a localized  area  in  the  interior  of  the  cylinder  of  liquid,  away  from 
walls  and  boundaries.  Most  measurements  were  made  at  a frequency  of 
36  kHz.  The  threshold  was  determined  by  increasing  the  acoustic  pressure 
amplitude  a small  amount,  waiting  for  a fixed  time  (normally  one  minute) , 
increasing  the  amplitude  and  so  forth  until  an  audible  snap  was  heard. 

The  pressure  amplitude  was  then  decreased  to  zero  for  a few  minutes,  and 
the  measurement  repeated  several  times.  Fresh  samples  were  allowed  to 
"age"  before  measurement  and  data  that  had  a time  dependence  were 
discarded.  Measurements  of  the  physical  parameters  such  as  surface  tension, 
gas  content  and  temperature  were  measured  in  situ  by  a du  Nuoy  ring 
tensiometer,  an  oxygen  electrode  and  a bead  thermister,  respectively. 

The  acoustic  pressure  amplitude  was  measured  at  the  center  of  the 
cavitation  zone  by  a calibrated  hydrophone  as  a function  of  input  current 
to  the  transducer  and  the  cavitation  chamber  calibrated  in  terms  of  this 
input  current. 

In  order  to  test  the  theory,  measurements  were  made  of  the  variation 
of  the  incipient  cavitation  threshold  as  a function  of  dissolved  gas 
content,  temperature,  and  surface  tension.  These  data  are  shown  in 
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Figs.  3,  4 and  5,  respectively.  Equation  13  contains  three  constants 
that  can  not  be  specified  initially.  These  constants  are:  y,  a constant 
that  measures  the  ability  or  inability  of  diffusion  to  keep  pace  with  the 

rapidly  expanding  cavity;  6,  a measure  of  the  advancing  contact  angle; 

■) 

and  (j>,  a measure  of  the  hysteresis  between  the  advancing  and  receding 
contact  angles.  Equation  13  was  then  applied  to  the  data  in  Figs.  3,  4 
and  5 and  a best  fit  obtained.  The  values  of  the  constants  obtained  were 
y = 1.0,  5 = 0.038  and  4*  = 35°.  The  implications  of  these  constants  will 
be  discussed  later. 

As  a further  result,  ordinary  tap  water  was  filtered  with  a 
Millipore  filter  with  a pore  size  of  0.45y  and  the  particulate  matter 
examined  with  a scanning  electron  microscope.  Figures  6 and  7 show 
photomicrographs  of  the  particulate  matter  examined.  Note  the  tremendous 
number  of  sites  that  are  available  for  gas  entrainment — one  particle 
such  as  this  could  provide  nuclei  for  thousands  of  cavitation  events. 

It  is  also  to  be  noted  that  the  crevices  are  small  angled,  in  support 
of  our  initial  hypothesis  of  cos  g •>  1. 

DISCUSSION 

We  have  produced  a theoretical  prediction  of  the  incipient 
threshold  for  acoustic  cavitation  inception  and  have  compared  our 
predictions  with  experimental  measurements.  In  Fig.  3 it  is  seen  that 
good  agreement  is  obtained  between  theory  and  experiment  for  the  variation 
of  the  cavitation  threshold  with  dissolved  gas  content.  We  have  found 
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that  a best  fit  required  y = 1.0;  an  interesting  result  in  that  this 
indicates  that  diffusion  of  gas  into  the  expanding  nucleus  is  rapid 
enough  to  keep  pace  with  the  expansion.  Perhaps  it  should  be  said  that 
the  expansion  is  slow  enough  for  diffusion  to  keep  pace.  This  agreement 
is  a welcome  result  for  several  observers  have  noted  the  dependence  of 
the  cavitation  threshold  on  the  gas  content  but  a correct  theoretical 
prediction  of  both  the  magnitude  and  the  dependence  of  the  effect  has  not 
been  previously  demonstrated. 

Figure  4 shows  the  dependence  of  the  cavitation  threshold  on  the 
surface  tension  of  the  liquid.  Most  cavitation  researchers13  would 
expect  a direct  dependence  of  the  threshold  on  the  surface  tension  due 
to  the  oft-used  term  2a/R  to  measure  the  pressure  in  the  interior  of  a 
cavity  of  radius  R due  to  a surface  tension  o . Hoyt8  has  even  suggested 
that  one  would  expect  the  analysis  of  Apfel4  to  predict  such  a direct 
dependence.  The  experimental  result  is  promising  for  it  suggests  that 
cavitation  reduction  may  result  from  addition  of  surfactants  to  the 
liquid  to  be  cavitated.  It  is  seen  that  an  increase  in  the  threshold  by 
a factor  of  nearly  four  is  experimentally  realizable.  The  author  plans 
to  suggest  this  effect  as  an  explanation  for  the  reduction  in  the 
cavitation  index  observed  by  Hoyt,8  and  others,  in  liquids  containing 
drag-reduction  agents,  which  reduce  the  surface  tension  considerably. 

It  is  seen  in  Fig.  4 that  the  inverse  dependence  is  predicted  by  the  theory 
as  well  as  the  correct  magnitude  of  the  threshold  for  two  values  of  the 
equilibrium  gas  pressure. 
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A further  test  of  the  applicability  of  the  theory  is  its  ability 

to  describe  the  variation  with  temperature.  Figure  5 shows  the  variation 

with  temperature  obtained  for  a frequency  of  22  kHz  and  an  equilibrium 

gas  concentration  of  0.35  bar.  Also  shown  are  the  normalized  results  of 

I 

Galloway14  who  had  calibration  problems  and  lists  values  of  the  threshold 
at  least  an  order  of  magnitude  higher  than  others.  His  measured 
dependence  of  the  threshold  with  temperature  is  seen  to  be  confirmed 
by  the  theory  as  well  as  our  own  measured  variation.  It  is  important 
to  note  that  in  order  to  obtain  the  correct  theoretical  dependence,  both 
the  vapor  pressure  and  the  equilibrium  gas  pressure  must  be  expressed  as 
a function  of  temperature.  We  have  used  the  following  variations  with 
temperature 

Pv  =(1.61  -(1.03  x 10-1)T  + (6.25  x 10-3)T2)xl04  dyn/cm2 

and 

Pg  = P“  (0.50  + 0.02T), 

where  Pg  is  the  equilibrium  gas  concentration  at  25°C  and  T is  in  C. 

It  should  be  noted  that  the  major  effect  of  temperature  is  through 
Pg  rather  than  Pv. 

It  is  noted  in  Eq.  13  that  the  threshold  is  independent  of  the 
frequency.  For  the  kilohertz  range  of  frequencies,  there  seems  to  be 
little  dependence  on  the  frequency2  and  thus  our  result  predicts  the 
correct  frequency  dependence. 
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igure  5-  Variation  of  the  cavitation  threshold  of  water  with  ter.perature  for 
two  values  of  the  equilibrium  gas  pressure.  The  curves  are  for  Eq.  13 
with  =1.0,  S =0.038  and  fi  =35°»  The  frequency  was  22  kHz  and  the  surface 
tension  was  70  dyn/cm.  The  plus  signs  are  the  normalized  values  of 
GallowaylH;  the  circles  are  measurements  from  this  study. 


In  order  to  obtain  the  correct  magnitudes  of  the  theory  to  agree  with 
the  results  of  Figs.  3,  4 and  5 it  was  noted  that  the  value  of  the  constants 
required  were  y = 1.0,  6 = 0.038  and  <J>  = 35°.  It  was  noted  earlier  that 

Y = 1.0  indicated  that  no  diffusion  correction  was  required.  6 = 0.038 

) 

implies  that  cos  (aA  - B)  = 0.038.  Thus  aA  = B + 88°,  and  if  we  assume 
that  B ~ 15°,  then  aA  *•  103°.  Bargeman  and  Van  Voorst  Vader12  have 
obtained  a value  of  106.5°  for  the  advancing  contact  angle  of  distilled 
water  on  parafin,  in  excellent  agreement  with  our  prediction.  The  other 
constant  used  in  the  analysis  is  <f>  which  is  given  by  <J>  = ajj  + B where 
01^  is  the  degree  of  hysteresis  between  the  advancing  and  receding 
contact  angle.  Again,  if  we  allow  B 15°,  then  a^  * 20°.  Furmidge11 
has  measured  the  angle  of  hysteresis  for  surfactant  plus  water  on  parafin, 
beeswax  and  cellulose  acetate.  His  measurements  give  a mean  hysteresis 
angle  of  26.7°,  which  agrees  with  our  required  angle  of  20°. 

We  can  use  these  constants  also  to  obtain  an  estimate  of  the  size  of 
nuclei.  Equation  4 gives  the  half-width  of  the  crevice  mouth.  Using  the 
values  o = 70  dyn/cm,  Py  - Pv  - Pg  » 1 x 106  dyn/cm2  and  cos  (aA  - B)  = 

0.038  we  find  that  the  full  width  of  a nucleation  site  in  a crevice  would 
be  approximately  0.05^m.  The  scanning  electron  micrograph  in  Fig.  7 shows 
that  the  tops  of  the  crevices  and  cracks  on  a mote  are  on  the  order  of 
0.5 <m.  Thus,  nucleation  sites  would  arise  from  air  trapped  in  the  bottoms 
of  some  of  these  crevices. 

In  summary,  we  have  measured  the  variation  of  the  incipient  cavitation 
threshold  in  water  with  such  variables  as  dissolved  gas  content. 
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Figure  6,  Scanning  electron  photomicrograph  of  a mote.  The  particle  was  removed 
by  filtration  of  distilled  water  with  a 0.45  ym  pore  size  filter,  shown  as 
the  mottled  area  in  the  background.  The  solid  horizontal  bar  in  the  fore- 
ground is  a 10  ym  marker.  Note  the  ragged  appearance  of  the  surface  showing 
the  numerous  possible  nucleation  sites. 
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Figure  7.  Scanning  electron  photomicrograph  of  a mote.  Enlarged  view  of  the 
mote  shown  in  Figure  6 • The  horizontal  bar  in  the  foreground  is  4 ym. 


surface  tension  and  temperature.  We  have  modified  the  analysis  of  Apfel 
to  obtain  an  expression  that  correctly  predicts  the  experimental  depended 
cies.  Finally,  our  model  of  nucleation  from  small-angled  crevices  and 
cracks  *n  solid  particulate  matter  is  substantiated  by  scanning  electron 

•I 

micrographs  of  filtered  particulate  matter. 
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PART  II 


A 


■I 

AIR  BUBBLE  GROWTH  IN  INSONATED  ROOT  TIPS  OF 
VICIA  FABA 


22 


INTRODUCTION 


r 

The  exposure  of  the  tips  of  root  tips  to  sound  irradiation  in  the 
megahertz  frequency  range  can  lead  to  growth  retardation  or  death  of  the 
tip,  depression  of  the  mitotic  index  and  induction  of  chromosomal 
anomalies.15  The  threshold  for  damage  at  1.0  MHz  is  a function  of  the 
ultrasonic  intensity  and  irradiation  time16  and  is  in  the  range  0.4  - 
0.7  W/cm2  for  irradiation  times  as  long  as  2 - 3 hours.  If  the  acoustic 
intensity  is  increased  to  levels  on  the  order  of  10  W/cm2,  the  root  tips 
often  are  killed.17  It  is  believed  that  the  damage  induced  by  the  ultra- 
sonic irradiation  is  the  result  of  the  activation  of  gas  bubbles  that  are 
present  in  the  root  tip.18  Large  amounts  of  intercellular  gas  exist  in 
a growing  root  system;  Gershoy18  et  al.,  have  examined  the  effects  of 
sound  irradiation  on  the  activation  of  this  gas.  They  have  observed  that 
the  intercellular  gas  has  caused  significant  acoustic  streaming  at 
intensities  on  the  order  of  0.35  W/m2,  with  gradual  increase  until  at 
intensities  of  35  W/m2,  free  mixing  of  cytoplasm  and  vacuolar  contents 
occurred  with  rapid  streaming  involving  the  entire  cell. 

Since  gas  bubbles  that  are  activated  by  sound  fields  emit  white 
noise,  subharmonics  and  harmonics  of  the  driving  frequency,  Coakley19 
has  examined  acoustic  emissions  from  insonated  root  tips  to  investigate 
the  role  of  activated  gas  bubbles  in  root-tip  damage.  His  system  was 
very  sensitive  to  acoustic  emissions  from  the  root  tip  during  irradiation 
in  that  the  root  tip  was  placed  near  the  center  of  a ring  transducer. 
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He  found  that  subharmonic  and  white  noise  emissions  occurred  at 
approximately  the  same  intensity  and  had  an  intriguing  history.  Figure  8 
below  is  a sketch  of  a typical  subharmonic  emission  from  a root  tip  at 
50  W/ cm2 . 
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Figure  8.  Typical  subharmonic  envelope, 
b- onset  of  subharmonic,  c- sound  off. 


a-sound  on. 


It  is  noted  that  there  was  a significant  delay  time  before  the  onset  of 
the  subharmonic  signal.  The  subharmonic  then  increased  in  intensity  quite 
rapidly  and  grew  to  a maximum  intensity.  There  was  a gradual  reduction  in 
intensity  until  after  a few  minutes  only  sporadic  activity  was  observed. 
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The  relatively  long  delay  time  between  the  initial  irradiation  and  the 
onset  of  the  subharmonic  looked  very  much  as  if  this  could  be  the  result 
of  the  growth  of  bubbles  by  rectified  diffusion  to  near  resonance  size, 
where  they  would  pulsate  nonlinearly  and  emit  a significant  amount  of 
radiation.  The  variation  of  the  delay  time  with  intensity  seemed  to 
confirm  this  speculation  as  there  was  a significant  dependence  of  the 
delay  time  on  irradiation  intensity,  as  shown  in  Fig.  9.  An  attempt  was 
made  to  explain  this  delay  time  in  terms  of  the  growth  of  small  nuclei 
by  rectified  diffusion. 


RESULTS 


A.  Applicable  equations 

The  equation  governing  the  rate  of  change  of  the  radius  of  an  air 
bubble  with  time  is  given  by  Eller"10  to  be 


dR  _ Dd  f2(l  - B2/s)  A-*2  2a -i 
dt  R '-3(1  - 32)2  ' '•pJ  ~ P . R'  ’ 


where  D is  the  diffusion  constant,  d is  the  ratio  of  the  concentration, 
in  mass  per  unit  volume,  of  dissolved  gas  to  the  density  of  the  gas  in 
equilibrium  with  the  solution.  Pa  the  acoustic  pressure  amplitude,  P0 
the  ambient  pressure,  and  a the  surface  tension.  The  resonance  factor 
3 2 is  given  by 
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Figure  9*  Variation  of  the  delay  time  for  noise  emmission  with  acoustic 
intensity  for  irradiated  root  tips  (results  of  Coakley1^). 


B2  - pU)2R2/[3n(P0  - pv)  + (3n  - l)2o/R]  , 


(15) 


where  <±>  is  the  angular  frequency  and  q is  a constant  that  measures  the 
thermal  behavior  of  the  gas  in  the  interior  of  the  bubble.  For  isothermal 
behavior,  r|  = 1.0;  for  adiabatic  cases,  q = 1.4.  It  was  discovered  that 
for  air  bubbles  driven  near  or  below  the  resonance  size  at  1 MHz,  q is 
within  2%  of  being  1.0.  Since  it  was  desired  to  integrate  this  equation 
from  nuclei-type  sizes  to  the  resonance  size,  it  was  necessary  to  account 
for  behavior  near  resonance.  Equation  14  is  undefined  when  g2  = 1,  the 
resonance  condition.  This  equation  was  modified  by  adding  a damping  term 
that  prevents  the  equation  from  becoming  undefined  at  resonance.18 
Thus,  the  equation  that  was  used  to  calculate  growth  times  was 


dR  Dd  f2(l  - B2/b)  r!V  20) 

dt  = R '3[(1  - m2  + V]  lP0J  ' P0R;  ’ 


(16) 


where  (2  is  a complicated  function  of  thermodynamic  constants  and  is 
discussed  in  more  detail  elsewhere.21 


B.  Theoretical  results 

It  was  discovered  that  for  these  large  acoustic  intensities, 
bubbles  would  grow  from  very  small  sizes  to  resonance  size  within 
relatively  short  times.  By  computer  integration  of  Eq.  16  it  was  found 
that  the  time  required  to  grow  from  an  initial  radius  of  5 x 10-6  cm  to 
resonance  size,  which  is  about  3.2  x 10-4  cm,  was  on  the  order  of  1 ms 


for  the  lowest  intensity  found  to  cause  subharmonic  emission 
(approximately  20  W/cm2).  Conditions  within  the  root  tip  are  obviously 
different  from  the  ideal  conditions  specified  for  growth  by  rectified 
diffusion  but  four  orders  of  magnitude  seem  too  much  to  make  up.  It  was 
concluded  that  the  delay  time  observed  between  the  commencement  of  the 
irradiation  and  the  onset  of  the  subharmonic  was  very  likely  not  growth 
by  rectified  diffusion. 

It  was  discovered,  however,  that  rectified  diffusion  effects  appeared 
to  be  present.  Figure  10a  shows  a fast  trace  scan  of  the  subharmonic 
emission  output.  There  appeared  in  each  trace  several  growth  sequences 
that  had  relatively  consistent  growth  times.  In  this  particular  scan, 

10  growth  sequences  can  be  observed  with  an  average  growth  time  of  0.3  ms. 
These  times  are  on  the  right  order  of  magnitude  for  rectified  diffusion 
growth  rates  so  a closer  examination  of  the  problem  was  made. 

Figure  11  shows  computer  integration  curves  of  the  time  required 
to  grow  from  a range  of  initial  sizes  to  resonance  size  for  various 
values  of  the  acoustic  intensity.  An  interesting  aspect  of  the  curves  is 
a leveling  off  of  the  times  to  a maximum  value  for  initial  radii  less  than 
about  0.5pm.  This  asymptotic  behavior  is  of  general  interest,  and  can 
be  demonstrated  by  studying  the  differential  equation  directly. 

Equation  16  is  a complicated  equation  near  resonance  but  if  it  is  examined 
far  enough  away  from  resonance,  it  can  be  integrated  directly.  If  it  is 
assumed  that  8 2 « 1,  then  Eq.  16  can  be  written  as 
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Figure  10a.  Oscilloscope  trace  of  subharmonic  emissions  from  an  insonated  root  tip 
The  sweep  rate  was  0.5  ms /cm  and  the  acoustic  intensity  was  about  30  W/cm2 . 


Figures  10b  & 10c.  Oscilloscope  trace  of  subharmonic  emissions  from  an  insonated 

root  tip.  The  sweep  rate  was  0.5  ms/cm  and  the  acoustic  intensity  was  about 
30  W/cm2.  Figure  b was  taken  immediately  after  the  onset  of  subharmonic 
emission;  Figure  c was  taken  approximately  100  seconds  later. 
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INITIAL  RADIUS  OF  BUBBLE  (jxm) 

Figure  11.  Growth  times  required  to  grow  from  an  initial  size  to  resonance 
size  for  various  values  of  the  acoustic  intensity.  Curves  are  for 
free  bubbles  in  saturated  water. 


dR  a Dd_  (2  r 2>\  zupq 
dt  “ R 3 Pn  J " PaR2 


If  the  substitutions  A » ^(^)2  — , and  B 


are  made,  then  this 


equation  can  be  written  as 


rT,  fR  R2dR 
* dt  = / AR=B  ’ 


which  can  be  integrated  by  parts  to  give  an  estimate  of  the  time  At 


required  to  grow  from  an  initial  radius  R0  to  a final  radius  R 


(that  is  significantly  less  than  Rres)  to  be 


^ - £*[4-  <r-r„>!  + <r-r„>  + t MSEiMi 


B o ..  fB  ARq 


A l2B 


This  equation  shows  two  interesting  aspects.  First,  for  certain  values 


of  the  initial  radius,  viz.,  Rt  -5.  B/A,  an  infinite  time  is  required  for 


2 r^Ai 2 2a 

growth.  But  Rt  = B/A  implies  that  -r-f— J = ~p  ~D~  and  this  is  essentially 

B *oKt 


the  threshold  for  rectified  diffusion.  For  intensities  lower  than  the 


threshold  value,  the  bubble  will  dissolve.  For  our  lowest  intensity  of 


interest,  I = 20  W/cm2,  Rt  = 3.5  x 10-6cm.  Note  that  slightly  above  this 


value  the  second  and  third  terms  in  the  equation  are  relatively  unimportant 


due  to  the  small  value  of  B/A  (*■  10“6)  and  the  time  is  essentially 


governed  by  the  simple  relationship 


..  - “a-1'  - ajifs*  «•>”.>■ 


Thus,  Che  asymptotic  approach  of  At  to  a constant  value  for  low  R is 
easily  shown  and  is  given  by 


At 


3Rp 2 rPjn  2 
4Dd  '•PA;  ’ 


(21) 


with  the  caveat  that  R > 3aP0/P^2.  This  result  asserts  that  a relatively 
large  range  of  bubble  sizes  will  grow  to  resonance  size  in  the  same  time 
interval,  an  interesting  and  important  deduction. 

It  is  also  of  interest  to  examine  the  behavior  of  air  bubbles  that 
may  be  larger  than  resonance  size.  Consequently,  Eq.  16  was  numerically 
integrated  to  show  the  time  required  for  a bubble  to  grow  from  an  initial 
radius  to  its  maximum  size.  Figure  12  shows  the  growth  of  a bubble  with 
initial  size  R = 5.0  x 10“6.  Note  that  the  bubble  grows  progressively 
more  rapidly  as  it  approaches  resonance  but  then  grows  steadily  until 
it  approaches  asymptotically  a maximum  size.  If  Eq.  16  is  examined 
directly,  it  is  seen  that  if  the  surface  tension  term  is  neglected,  then 
the  growth  rate  becomes  zero  for  B2  = 8,  and  negative  for  larger  values 
of  B2.  Since  B2  = 1 is  the  resonance  condition,  it  is  seen  that  air 
bubbles  with  radii  larger  than  /T  times  their  resonance  size  will 
actually  be  forced  to  smaller  sizes  by  rectified  diffusion. 

In  Fig.  10b  and  c are  two  oscilloscope  traces  of  the  emitted 
subharmonic  output  from  an  insonated  root  tip.  The  first  trace  was  taken 
immediately  after  the  acoustic  emission  commenced;  the  second  was  taken 
nearly  100  seconds  later.  The  first  trace  shows  a few  relatively  long 
and  relatively  consistent  growth  times;  the  second  shows  several  relatively 
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Figure  12.  Life  history  of  growth  of  bubble.  This  is  for  an  intensity 
of  30  W/cm^  and  a frequency  of  1 MHz. 


short  and  relatively  inconsistent  growth  times.  These  traces  suggest,  in 
the  light  of  the  above  theoretical  analysis,  that  air  bubbles  may  be  growing 
to  resonance  size  from  nuclei,  and  then  breaking  up  into  a group  of  smaller 
bubbles  by  surface  wave  activity. 


DISCUSSION 


The  analysis  of  acoustic  emissions  from  insonated  root  tips  indicated 
that  it  is  unlikely  that  the  relatively  long  delay  time  of  seconds 
between  the  commencement  of  the  irradiation  and  the  onset  of  subharmonic 
and  white  noise  is  due  to  growth  of  air  bubbles  by  rectified  diffusion. 
Theoretical  examination  of  growth  from  a range  of  nuclei  sizes  indicates 
growth  times  on  the  order  of  tenths  of  milliseconds.  It  is  possible  that 
certain  constraints  may  exist  within  the  root  system  that  would  slow  down 
rectified  diffusion  but  none  seem  likely  to  affect  it  to  four  orders  of 
magnitude.  Further,  since  oscilloscope  scans  of  the  subharmonic  emissions 
show  growth  times  of  the  radiated  signal  on  the  order  of  0.3  ms,  and  since 
calculated  growth  times  for  rectified  diffusion  are  in  the  range 
0 to  1.3  ms,  it  seems  likely  that  air  bubbles  are  growing  by  rectified 
diffusion  but  at  the  shorter  times  rather  than  the  longer.  Unstable 
growth  of  cavities  through  vaporous  cavitation  can  occur  in  these  growth 
times  also  and  must  be  considered  as  an  alternative  explanation.  Its 
examination  is  too  complicated  to  consider  here  however. 

There  is  a mechanism  that  will  lead  to  delay  times  on  the  order  of 
10  seconds  and  to  subsequent  violent  activity.  The  time  required  for 
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air  bubbles  much  larger  than  resonance  size  to  reduce  in  size,  is  on  the 


1 

order  of  seconds.  For  example,  for  an  intensity  of  30  W/cm,  three  seconds 
is  required  for  an  air  bubble  of  10u  radius  to  reduce  to  a radius  of 
7.9pm,  the  radius  for  which  rectified  diffusion  ceases  at  a frequency  of 
1 MHz;  the  time  required  for  the  radius  to  change  from  20pm  to  8.3pm  is 
17  sec.  Measurements  of  the  intercellular  spaces  are  on  this  order  of 
magnitude.16  Thus,  it  is  possible  that  this  intercellular  gas  is  reduced 
in  size  until  it  can  develop  surface  waves.  These  surface  waves  cause 
the  bubble  to  break  up  into  several  small  bubbles,  which  grow  by  rectified 
diffusion  to  resonance  size  again  and  again. 
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800  No.  th  Quincy  Street 
Arlington,  Virginia  22217 

Naval  Research  Laboratory  3 copies 

Department  of  the  Navy 
Attn:  Technical  Library 
Washington,  D.  C.  20375 

Office  of  the  Director  of  Defense  3 copies 

Research  and  Engineering 
Information  Office  Library  Branch 
The  Pentagon 
Washington,  D.  C.  20301 

U.  S.  Army  Research  Office  2 copies 

Box  12211 

Research  Triangle  Park 
North  Carolina  27709 

Defense  Documentation  Center  12  copies 

Cameron  Station  (TC) 

Alexandria,  Virginia  22314 

Director,  National  Bureau  of  Standards  1 copy 

Attn;  Technical  Library 
Washington,  DC  20234 

Commanding  Officer  3 copies 

Office  of  Naval  Research  Branch  Office 
538  South  Clark  Street 
Chicago,  Illinois  60605 
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3 copies 


Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91101 

San  Fra  .cisco  Area  Office  3 copies 

Office  of  Naval  Research 
One  Hal  1 Idle  Plaza 
Suite  601 

San  Francisco,  California  94102 

Commanding  Officer  3 copies 

Office  of  Naval  Research  Branch  Office 

666  Summer  Street 

Boston,  Massachusetts  02210 

New  York  Area  Office  1 copy 

Office  of  Naval  Research 
715  Broadway,  5th  Floor 
New  York,  New  York  10003 

Director  1 copy 

U.  S.  Army  Engineering  Research 
and  Development  Laboratories 
Attn:  Technical  Documents  Center 
Fort  Belvoir,  Virginia  22060 

0DDR6E  Advisory  Group  on  Electron  Devices  3 copies 

201  Varick  Street 

New  York,  New  York  10014 

Air  Force  Office  of  Scientific  Research  1 copy 

Department  of  the  Air  Force 
Bolling  AFB,  D.  C.  22209 

Air  Force  Weapons  Laboratory  1 copy 

Technical  Library 
Kfrtland  Air  Force  Base 
Albuquerque,  New  Mexico  87117 

Air  Force  Avionics  Laboratory  1 COpV 

Air  Force  Systems  Command 

Technical  Library 

Wr ight-Patterson  Air  Force  Base 

Dayton,  Ohio  45433 

Lawrence  Livermore  Laboratory  • 1 copy 

Attn:  Dr.  W.  F.  Krupke 
University  of  California 
P.  0.  Box  808 

Livermore,  California  94550 
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f 


Harry  Diamond  Laboratories 

Technical  Library 

2800  Powder  Mill  Road 

Adelphl,  Maryland  20783 

1 copy 

Naval  Air  Development  Center 

Attn:  Technical  Library 

Johnsvi 1 le 

Warminster,  Pennsylvania  18974 

1 copy 

Naval  Weapons  Center 

Technical  Library  (Code  753) 

China  Lake,  California  93555 

1 copy 

Naval  Training  Equipment  Center 

Technical  Library 

Orlando,  Florida  32813 

1 copy 

Naval  Underwater  Systems  Center 

Technical  Library 

New  London,  Connecticut  06320 

1 copy 

Commandant  of  the  Marine  Corps 

Scientific  Advisor  (Code  RD-l) 

Washington,  DC  20380 

1 copy 

Naval  Ordnance  Station 

Technical  Library 

Indian  Head,  Maryland  20640 

1 copy 

Naval  Postgraduate  School 

Technical  Library  (Code  0212) 

Monterey,  California  93940 

1 copy 

Naval  Missile  Center 

Technical  Library  (Code  5632.2) 

Point  Mugu,  California  93010 

1 copy 

Naval  Ordnance  Station 

Technical  Library 

Louisville,  Kentucky  4021 A 

1 copy 

Commanding  Officer 

Naval  Ocean  Research  6 Development  Activity 
Technical  Library 

NSTL  Station,  Mississippi  39529 

1 copy 

Naval  Explosive  Ordnance  Disposal  Facility 
Technical  Library 

Indian  Head,  Maryland  20640 

1 copy 

AO 


L 


Naval  Ocean  Systems  Center 

Technical  Library 

San  Diego,  California  92152 

1 

copy 

Naval  Surface  Weapons  Center 

Technical  Library 

Dahlgren,  Virginia  22448 

1 

copy 

Naval  Surface  Weapons  Center  (White  Oak) 

Technical  Library 

Silver  Spring,  Maryland  20910 

1 

copy 

Naval  Ship  Research  and  Development  Center 

Central  Library  (Code  L42  and  L43) 

Bethesda,  Maryland  20084 

1 

copy 

Naval  Avionics  Facility 

Technical  Library 

Indianapolis,  Indiana  46218  

1 

copy 
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